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Abstract 



We consider global geometric properties of a codimension one manifold embedded 
in Euclidean space, as it evolves under an isotropic and volume preserving Brownian 
flow of diffeomorphisms. In particular, we obtain expressions describing the expected 
rate of growth of the Lipschitz-Killing curvatures, or intrinsic volumes, of the manifold 
under the flow. 

These results shed new light on some of the intriguing growth properties of flows 
from a global perspective, rather than the local perspective, on which there is a much 
larger literature. 



1 Introduction 

We are interested in Brownian flows $ st , < s < t < oo from R" — > ]R", obtained by solving 
the collection of stochastic differential equations 



where we write $t for $ot when there is no danger of confusion. Here, d denotes the 
Stratonovich stochastic differential and Ut (x) is a vector field valued Brownian motion with 
smooth spatial covariance structure, on which wc shall have more to say in the subsequent 
section. However, we note already that we shall assume U is such that, with probability 
one, for each s < t, 

(i) & s t is a C 2 diffeomorphism. 
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(ii) $ s t is volume preserving; i.e. for any compact D C R™, A„($t(D)) = X n (D), where 
A„ is Lebesgue measure in R™. 

(iii) $ s t is isotropic in the sense of (8) below. 

It is standard fare, following from (1) and our three assumptions, that $ ut o $ stl = $ st , that 
$tt is the identity map on R™, that Q s t(x) and <&~ t (x) are j° m tly continuous in x, s, t as are 
the spatial derivatives 

and C$7* 1 C 33 )) an d that the 'increments' ^ Sl t! , <&s 2 t 2 j ■ ■ ■ > ^s^tn are independent for all si < 
ii < S2 < ^2 < • • • < s„ < t n . For a full study of isotropic Brownian flows, with history and 
references, we refer the reader to Kunita's monograph, [6]. 

The study of the evolution of curvature under such flows was pioneered by Le Jan in [8] , where 
he established the positive recurrence of the curvature of a curve moving under an isotropic 
Brownian flow. Quite recently, Cranston and LeJan [2] followed this with a striking analysis 
of the growth of local curvature. Working with isotropic and volume preserving flows in R", 
they took a codimcnsion one manifold AI embedded in R™, and considered its image under 

the flow, which we denote by M t = $ t (M). 

Taking a point x S M, they developed an Ito formula for the symmetric polynomials of 
the principal curvatures of <!>t(Af) at the points $> t (x), including the mean and Gaussian 
curvatures. In addition, and this will be more important for us, they showed that these 
polynomials grow exponentially in time, with mean rates that are related to the Lyapunov 
exponents of the flow. In simple terms, this means that the manifold M tl while it may begin 
at time zero as something as simple as the unit sphere S 1 ™ -1 in R n (which has unit Gaussian 
curvature everywhere) it tends to develop sharply rounded 'corners' as time progresses. 

A somewhat different set of results can be found in a series of papers [3, 4] authored by 
Cranston, Scheutzow and Steinsaltz. In particular, the combined results of [3, 4, 11] show 
that, for an isotropic Brownian flow with n > 2, there arc positive constants c and C such 
that for each compact and connected set D C R™ with at least two points, 

c < inf sup liminf -(<£> t (x) 7 u) < limsup - sup ||$t(x)|| < C, (3) 

ueS"- 1 x eD t^>oo t t^oo t x£D 

almost surely if the top Lyapunov exponent is strictly positive, and with strictly positive 
probability, otherwise. 

One implication of this result is that while $ t (M) is homotopically equivalent to M, for 
large t it will 'look' quite different. One way to measure this difference, at a global level, 
is via their Lipschitz-Killing curvatures. Since M has dimension (n — 1) there are n such 
curvatures, Co(M t ), £i(M t ), . . . , C n -\(M t ). The first of these, Co(M t ), is the Euler-Poincare 
characteristic of M t , which, because M and M t are homotopically equivalent, is the same 
as that of M, and so independent of t. The last of these, £„_i(M t ), gives the (n — 1) 
dimensional surface measure of M t and most definitely docs change with time, as do all 
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the remaining Cj(M t ). Further information on the geometric roles of the Lipschitz-Killing 
curvatures is given in the following section. 

In view of the results of Cranston, Scheutzow and Steinsaltz described above, one would 
expect that the Cj(M t ), 1 < j < n — 1, would grow rapidly in time, as parts of the set M t 
begin to stretch in various directions at rate t. That this is indeed the case is a consequence 
of the following theorem, one of the two main results of this paper. 

Theorem 1.1 Let AI be a smooth codimension one manifold embedded in R n and M t its 
image under $ t; where $ s t is an isotropic and volume preserving Brownian flow of C 2 
diffeomorphsims o/R™. Then, forO < k < n — 1, the expected rate of growth of the Lipschitz- 
Killing curvatures is given by 

/ 1 1 „ . . f(n- k - l)(n+ l)(k + l)/x 2 A 
E (M t )} = £ n _ fe _ 1 (M)ex P n £ + ^ ^-J , (4) 

where [i2 is the second moment of the spectral measure F of (5). 

The proof of Theorem 1.1 relies on the fact that, loosely speaking, for 1 < k < (n — 1), the 
(n — k — l)-th Lipschitz-Killing curvature of a manifold can be obtained as an average, over 
the manifold, of the fc-th order symmetric polynomial of the principal curvatures. However, 
as we have already noted above, these have been studied in detail by Cranston and LeJan 
[2]. Consequently, our proof relies very heavily on their paper, to the extent that one could 
consider this paper as an addendum to theirs. Nevertheless, we believe that the results 
are of independent interest, in that they lift the local approach of [2] to a global scenario. 
Finally the case k = 0, which corresponds to the (n — 1) dimensional surface measure of the 
manifold or cquivalently the (n — l)-th Lipschitz-Killing curvature, is a simple consequence 
of Lemma 4.1. 

The remainder of this paper is organized as follows. In Section 3 we provide the precise 
definition of Lipschitz-Killing curvatures and the required geometric background, followed 
by the proofs of the main results of the paper in Section 4. 



2 Brownian flows 

This section is not so much about Brownian flows per se, for which we refer you back to 
the references of the Introduction, but rather about setting up notation. Since we plan to 
use the main result of [2] to prove Theorem 1.1, we shall adopt the notation of that paper 
without much explanation. You can find missing explanations in [2]. 

The first step is to define the vector field valued Brownian motion U driving the flow in 
(1). We take this to be a zero mean Gaussian process from R+ x M" to R" with covariance 
structure given by 

E{U?(x)U l s (y)} = {tAs)C kl (x-y), l<k,l<n, 
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where each C can be written in the form 

C kl (z) = / / e^^-AVn-i^W, (5) 
Jo JS"- 1 

for a normalized Lebesgue measure cr„_i on S™ -1 and a non-negative measure F on R + , 
with m-th moment, /it m , given by /i m = /„ F(dp). 

The various spatial derivatives of U , implicitly assumed to exist, are denoted by 

W s = — — , B) 



Writing (■, •} for quadratic covariation, it is not hard to check that 

{dWj(t, y ),dW?{t, y )) = — ^[(n + l)4^-^f-^]di, (6) 
{dB} k {t,y\dWV{t,y)) = 0, (7) 

for any 1 < k, l,p,q< n, implying E(^2 Wl) 2 = 0, hence volume preserving, and 

((dB(u, u), v), (dB(u, u), v)) = n(n+ 3 2 ^ n + 4) [(n + 3) jM| 4 |M| 2 - 4{u, vf\\u\\ 2 ] dt, 
for all vectors u, v G R". 

This particular choice of the covariance function makes the flow isotropic in the sense that 
the spatial covariance matrices C(x) = (C fe '(x))™ fe=1 satisfy 

C(x) = G*C(Gx)G (8) 

for any real orthonormal matrix G, as well as making the flow volume preserving. 

With the flow defined, we now turn to setting up the notation required for studying its 
(differential) geometry. Our basic references for this are Lee [9, 10] and Part II of [1]. 

We start with a codimension one Riemannian manifold M embedded in R" and, for x G M, 
let u = {ui}^^ be an orthonormal basis of T X M, the tangent space at x G M. (Note 
that u actually depends on x, but we shall not write this explicitly.) Then by a simple 
push-forward argument, D$ t {x)ui = Uj(t) G T Xt M t , where xt = &t{x) (cf. (1)) and D$t(x) 
is as defined in (2). Furthermore, 

dui(t) = dWui{t) = dWui(t). 



Writing lit : T Xt M. n — > T Xt M t as the orthogonal projection onto T Xt M t , and V for the 
canonical connection on the ambient Euclidean space R™, the second fundamental form at 
Xt G M t is given by 

S t (u(t),v(t)) = (I-U t )V u(t) v(t), 
for any u(t),v(t) G T Xt M t . Subsequently, scalar second fundamental form is defined as 

S„ t (u(t),v(t)) = (5 t («(t), «(*)), i/ t ), 
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where v t denotes the unit normal vector field and {•,•) now denotes the usual Euclidean 
inner product rather than quadratic covariation. 

It is standard fare that the scalar second fundamental form can be used to induce a linear 
operator on the exterior algebra K k (T Xt M t ) of alternating covariant tensors on T Xt M t , built 
over the usual wedge product, for each 1 < k < (n — 1). First we define Sj>t as 

k 

Sl k ;\u ll (t)A...Au lk (t), u mi (t)A...Au mk (t)) = ^(-l)""n 5 ^( M ^)W' u ^W)' 

o-es fc 3=1 

where {ui p (t)}, {u mq (t)} C {ui(t)}\Z~i 1 \ &k is the collection of all permutations of {1, . . . , k} 
and I],? denotes the sign of the permutation a. This gives rise to the linear operator 

u ll (t)A...Aui k (f) « SW(u h (t) A...Aui fc (t),-), 

where {u,(t)} |"7 1} C T Xt M t is a basis of T Xt M t . 

The last and the most important remaining definition is that of the trace of St£ . For this, 
however, we need some more notation. For 1 < k < (n — 1) define the index set Ik by 

h = {m e {1, . . . , n - l} k : mi < mi < ■ ■ ■ < mk} ■ 
Then, for I € Ik, define 

\T\ = h + --- + i k , 

aft) = ui 1 (t) A • • • A ui k (t), 

a\t) = {-l) W+k Ul (t) A---Au h (t) A--- Aui k (t) A • • • A 
a(t) = Ui(t) A • • • A u n -i(t), 

where the hatted vectors are understood to be omitted from the wedge product. Now, for 
l,rn £ Ik, define 

(afi),<m(t)) = det((u k (t),u mj (t)}), (9) 
and naturally \\a(t) || 2 = det ({ui(t), u 3 (t))). 

We now have all that we need to define the all important trace, TrS 1 ^ , as 

TrSW = {al ^ { Q t)} S( k Haft),a A m (10) 

where the Einstein summation convention is carried over the indices l,m € Ik- Now for the 
case k = 0, which has thus far remained untouched, we define TrS^ = 1. 

The temporal development of this trace was studied in detail in [2], where Cranston and 
LeJan proved that {TrS 1 ^}^ 1 ^ is a (n — l)-dimensional diffusion. This is the precise 
formulation of the result we were referring to in the Introduction when we spoke of the Ito 
formula of symmetric polynomials of principal curvatures, which are essentially equivalent 
to the traces {TrS'^ ) }^™" 1 1) . 
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3 Lipschitz-Killing curvatures 



There are a number of ways to define Lipschitz-Killing curvatures, but perhaps the easiest 
is via their appearance in the so-called tube formulae, which, in their original form, are due 
to Weyl [12]. (For more details and applications see either the monograph of Gray [5] or 
Chapter 10 of [1].) 

To state the tube formula, let M be a C 2 , (n — l)-dimensional manifold embedded in M™ 
and endowed with the canonical Riemannian structure on R™ . The tube of radius p around 
M is defined as 

Tube(M, p) = {x e M" : d(x, M) < p} , 

where 

d{x,M) = inf 

Weyl's tube formula states that there exists a p c > 0, known as the critical radius of M, 
such that, for p < p c , the volume of the tube is given by 

n-l 

A„(Tube(M, j0 )) = Y,P n ~^n-jCj{M), (11) 

3=0 

where ujj is the volume of the j-dimensional unit ball and Cj (M) is the j'^-Lipschitz-Killing 
curvature of M. 

Writing Tij for j-dimcnsional Hausdorff measure, it is easy to check from (11) that C n -i (M) = 
Tln-i(M). That is, it is the surface 'area' of M. Cq{M) is the Euler-Poincare char- 
acteristic of M, and while the remaining Lipschitz-Killing curvatures have less transpar- 
ent interpretations, it is easy to see that they satisfy simple scaling relationships, in that 
Cj(aM) = a- ' £j(M) for all 1 < j < n — 1, where aM = {x G M" x = ay for some y £ M}. 
Furthermore, despite the fact that defining the Cj via (11) involves the embedding of M in 
R™, the Cj(M) are actually intrinsic, and so independent of the embedding space. 

While (11) characterizes the Cj(M) it does not generally help one compute them. There are 
a number of ways in which to do this, but we choose the following, which is most appropriate 
for our purposes, (cf. [1, 5] for further details and examples) 

Cn-k-xiM) = K n , k f f TrSWl^Mi-^Hoid^Hn-xidx), (12) 
Jm Js(r) 

where K n ,k = 2 ( 7r )( fc + 1 )/ 2 fc! ^ (2 / anc ^ N X M is the normal cone to M at the point x. Since 
M has codimension 1 in R™, each N X M contains only the outward normals to T X M and is 
of unit dimension. 

In general, we shall write S'(R n ) for the unit sphere in R™. Thus the S(M.) appearing in 
(12) contains only the two vectors +1 and —1 in R and Ho is counting measure, which 
makes the integral over S(M) a rather pretentious way of writing things, as, indeed, was the 
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introduction of the normal cone. Nevertheless, both will be of use to us later on, when we 
discuss possible generalisations of our results. 



We, of course, are interested in the temporal evolution of the Cj{M t ) and it follows directly 
from (12) that, with the notation of the previous section, 



C n - k -i(M t ) = K n . k / / TrSW\ N ^ Mt {- Vt )na(dvt)Hn-i(dxt) 



= K„ 



TrS%> yj det ( (u t (t ) , u 3 (t ) } ) 1 Nlt M t ( - v t ) H {dv t )H n -i (dx) 



= K n . k / / TrS^\\at\\lN Tt M t (-^)no(dut)n n -i(dx). (13) 

J M JS(K) 

4 An Ito formula for C n -k-i(Mt) 

Before commencing a serious stochastic analysis of (13) we recall some of the results and 
further notation from LeJan [7]. 

Let f(t) =6(«)A...A &(t) and xjj(t) = ip±(t) A . . . A rfi k (t), where {&(*)}, {^(t)} C T Xt M t . 
Then, by Lemma 3 of [7], 

dm,m) = E(Wf( f )^(*)> + m,r^(t)))dw;(t) + k{n ~ kh2 m,m)dt, 

L < J n 

1,3 

where 

rfm = e j A ^2(-l) l+1 Ut),e l )^(t) A • • • Afc(t) A • • • A&(t) j , 

with {e fc }£ =1 being the standard basis of ffi™, and ifVKO is defined similarly. 

It follows immediately from the above that if = ?/>(£) = a(i) = u\(f) A ■ ■ ■ A u„_i(i), 
where itj(i) = D^> t {x)ui and (iti, . . . ,it n _i) is an orthonormal basis of T X M, then 

n— 1 / _ I \ 

d||a(i)|| 2 = ||a(i)|| 2 (2^^(t) + ^— 

Now we derive an Ito formula for ||a(f)|| and use it together with Theorem A. 2 of [2] to 
obtain an expression for the Ito derivative of the Lipschitz-Killing curvatures. 

Lemma 4.1 Let M be a smooth (n — 1)- dimensional manifold embedded in W 1 and M t its 
image at time t under the stochastic, isotropic, and volume preserving flow $t described in 
Section 2. Then, in the notation of Section 2, 

d\\a(t)\\ = \\a(t)\\(XdWm + ( "~ n ^ n + + i )M2 ^ 

i— 1 
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Proof: Using the standard Ito formula and (6) we obtain 
d||a(t)|| = d(\\a(t)\\^ 

n-l 



i=l 
n-l 



-i)(n + r 

2n(n + 2) 



■wii(£^w + (n ^; ) "^ 



= 53[s(*- 1 )(a I ;(t),a^(t))(dB(«, p (t),« TO( (t)),^) 

n-l 

+ TrS« (fcdW£(t) - 2 ^ dW7(t) 



ll«(*)H s 



i=l 



□ 



We need just a little more preparation before we can turn to our main result. 
Let aft) = tiji A • • • A ui k (t), be a A:-form for Z S ifc, then define 

a fp {t) = (-l) p+1 u h {t) A-- - Ai, r (i) A-AU| t (i), (14) 

for 1 < A < (n — 1), where Z £ Ifc, Z p e /fc-i and 1 < p < k. 
Rewriting the above expression as 

«,;(*) = (-rMf'WA-Aui^w, (15) 

defines u[ . 

Then according to Theorem A. 2 of [2], for 1 < k < n — 1 

(c/(*),a™(i)) 



(n + l)fc(n — Zc)^2 



K*)ll s 

TrS«di, (16) 



2n(n + 2) 

where the Einstein summation convention is carried over the indices Z and to 
We now have everything we need to present the main result of this paper. 

Theorem 4.1 Retain the assumptions and notation of Lemma 4-1- Let Ck be the Lipschitz- 
Killing curvatures defined by (13). Then the Ltd derivatives of the Lip schitz- Killing curva- 
tures for 1 < k < n — 1 are given by 
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dC n -k-i(M t ) 



K. 



M J S 



( E Si k - 1 \a Cv {t),a^M{dB{u lp {t),u m M^t) 
x (o , (t),a' B (t)>||a(t)||- 1 

ra-1 

+ Tt5W]|a(t)||(fe^-^^) 

+ £ 5fW(a p a^)((^ r (i), a*(i)> + (afy), ^ a™(t)))dW^ ||a(i) r 1 



x lN xt M t {-vt)Ho(dvt)TLn-i{dx) 
, (n-fc-l)(n + !)(£; + l)/i 2 



2n(n + 2) 



C n -k-l(M t )dt, 



(17) 



where TLk is k- dimensional Hausdorff measure and N Xt M t is the normal cone to M t at 



Before giving the proof, we note that (17) simplifies considerably when k = 0, in which case, 
as we have already noted, C n -i{M t ) = H n -i(M t ). Then 

C n -i{M t ) = K(n,0) f I ||a(i)||liv St M t (-^)Wo(^t)W n _i(dz) 

JM JS(R) 

= f \\a(t)\\H n -x(dx). 

JM 

Hence, by Theorem 4.1, 

d£„_x(M t ) = f \\a(t)\\ (V dWt(t))n n ^(dx) + (n ~ + ( )/X2 £»_ 1 (M t )dt 

V fr^ 7 2n(n + 2) 

.EVi passant to the proof we remind the reader of the dependence of the various integrands 
appearing in (17) on the space parameter x £ M through the vector field U, its various 
spatial derivatives W and B, and the tangent vectors, at x £ M. 

Proof of Theorem 4.1: We start with 

d(TrSif\\a(t)\\) = (d(Tr5W)|| Q (t)||+Tr5W<||a(t)||) + (dTr5«,d||a(t)||) 
= I + 11 + III. 

We shall obtain a closed form expression for each of these terms. By (16) we have 
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(a T (t),a™{t))\Wt)\\- 1 



I = (d(TrSW)\\a(t)\\ 

k 

i.p—1 

71-1 

+ TrS« [MCW - 2 E dW?(t)] ||a(t)|| 
(n + l)k(n - k)n2 



2n(n + 2) 



-TrSW\\a(t)\\dt. 



(18) 



Using Theorem 4.1 we find 

II = TrS™d(\\a(t)\\) 



TrSW\\a(t)\\(j2dWm 



(n- l)(n + l)/i 2 
2n(n + 2) 



dt 



Finally using (6), (16) and Theorem 4.1 we have 
III 

[dTTS%>,d\\a{t)t 



(19) 



(20) 



= TrS« ||a(t)|| ( (kdW: -2^dWt),J2 dW t 

i-l i=l 

n-1 

+ £sW(aKt),o*(t))((7^^ 

(21) 



; = 1 



(n l)(k + 2)to^k> hmdt + 2(n-fe-l) M2 ^ )||aft)|| ^ 



n(n + 2) 



i(n + 2) 



Summing (18), (19) and (20) we have 
d(TrS«||a(t)||) 

= [^Sf 1) (a i ;(t) 1 ^(t))^(%(f),^(f)) 1 <' ( )](/(t) 1 a™(f))|ai|- 1 



+ TrSW\\a(t)\\ [kdWZ(t) - ]T dW?(t)' 

i=l 

+ X)^ ) (aKt),a*(«))((7f'a r (t),a fll (t)) + (afy.T? dW^(t) ||«(t) || - 1 

, (n-ft-l)(n + l)(Ar + l)Ais, 



2n(n + 2) 



■Tr^||a(t)||A. 
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Substituting the above in (13) gives (17) and so the theorem. 



□ 



We can now easily deduce Theorem 1.1 of the Introduction. 

Proof of Theorem 1.1: In (17), we note that with the single exception of the last term, all 
terms are zero mean martingales due to the presence of the martingale integrators dW {t) 
or dB- (t). Therefore, taking expectations in (17), after taking the integral over time t, 
immcdiatley yields 

E{L M (M t )} = + + f\ {Cn _ k _ l{ M s )} ds. 

2n[n + 2) J 

Solving this linear differential equation gives (4), and we are done. □ 

We close with one further result, and two open problems. 

As for the first open problem, we believe that, in the setting of Theorem 1.1, 

lim I lo ■ ( C n-k-i{M t ) \ (n-k-l)(n+l)(k + l) m 
t^ t ° 8 y C n -k-i(M) J 2n(n + 2) 

where the limit here is in L\. At this point we do not have an air tight proof of this. We 
would also like to add that the almost sure growth rates may well be different from the ones 
conjectured above. 

As for the additional result, recall that throughout the paper, we have assumed that M was 
a codimension one manifold in W l . From a technical point of view, this has a substantial 
simplifying effect on the definition (12) of Lipschitz-Killing curvatures. If dim(M) ~ m < 
(n — 1), then (12) changes in that the normal cones are now of dimension (n — m), 5(R) is 
replaced by S(M. n ~ m ) and so is of dimension (n — m—1), and Hq and 7i n -x are replaced by 
7i rl _ m _i and TL m , respectively. (The constant also changes, but this is less important. See 
[1] for details.) All told, we have 

£ m -k(M) = K' mjk [ f TrS?h NaM {-v)H n ^-x{dv)H m (dic), (22) 

J M JS(R™-™) 

for some constants K' m k . For all k ^ here, we have found that the complications introduced 
by increasing the dimension of the normal space are such that computations analogous to 
those we have carried out are forbiddingly complex. 

For k = the trace term in (22) disappears, and so it is not hard to show that 

d£ m (M t ) = I \\a(t)\\(f2dWi(t))H m (dx) + m(n ~ £ m (M t )dt. (23) 

Jm ' 2n(n + 2) 

Since C m (M t ) is the m-dimcnsional content of M t , this is a far from uninteresting result. 

Of course, as before, this implies that 

E { £ m (M t)} = £ ro (M )eX p ( ^^t) 1 ^ ) ■ ( 24 ) 
We have not, however, been able to find corresponding results for the more general case. 
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